Linking Spontaneous Behavioral Changes to Disease Transmission Dynamics: Behavior Change Includes Periodic Oscillation.

Behavior change significantly influences the transmission of diseases during outbreaks. To incorporate spontaneous preventive measures, we propose a model that integrates behavior change with disease transmission. The model represents behavior change through an imitation process, wherein players exclusively adopt the behavior associated with higher payoff. We find that relying solely on spontaneous behavior change is insufficient for eradicating the disease. The dynamics of behavior change are contingent on the basic reproduction number R a corresponding to the scenario where all players adopt non-pharmaceutical interventions (NPIs). When R a < 1 , partial adherence to NPIs remains consistently feasible. We can ensure that the disease stays at a low level or maintains minor fluctuations around a lower value by increasing sensitivity to perceived infection. In cases where oscillations occur, a further reduction in the maximum prevalence of infection over a cycle can be achieved by increasing the rate of behavior change. When R a > 1 , almost all players consistently adopt NPIs if they are highly sensitive to perceived infection. Further consideration of saturated recovery leads to saddle-node homoclinic and Bogdanov-Takens bifurcations, emphasizing the adverse impact of limited medical resources on controlling the scale of infection. Finally, we parameterize our model with COVID-19 data and Tokyo subway ridership, enabling us to illustrate the disease spread co-evolving with behavior change dynamics. We further demonstrate that an increase in sensitivity to perceived infection can accelerate the peak time and reduce the peak size of infection prevalence in the initial wave.

Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Estimating the population effectiveness of interventions against COVID-19 in France: A modelling study.

BACKGROUND: Non-pharmaceutical interventions (NPIs) and vaccines have been widely used to manage the COVID-19 pandemic. However, uncertainty persists regarding the effectiveness of these interventions due to data quality issues, methodological challenges, and differing contextual factors. Accurate estimation of their effects is crucial for future epidemic preparedness. METHODS: To address this, we developed a population-based mechanistic model that includes the impact of NPIs and vaccines on SARS-CoV-2 transmission and hospitalization rates. Our statistical approach estimated all parameters in one step, accurately propagating uncertainty. We fitted the model to comprehensive epidemiological data in France from March 2020 to October 2021. With the same model, we simulated scenarios of vaccine rollout. RESULTS: The first lockdown was the most effective, reducing transmission by 84 % (95 % confidence interval (CI) 83-85). Subsequent lockdowns had diminished effectiveness (reduction of 74 % (69-77) and 11 % (9-18), respectively). A 6 pm curfew was more effective than one at 8 pm (68 % (66-69) vs. 48 % (45-49) reduction), while school closures reduced transmission by 15 % (12-18). In a scenario without vaccines before November 2021, we predicted 159,000 or 168 % (95 % prediction interval (PI) 70-315) more deaths and 1,488,000 or 300 % (133-492) more hospitalizations. If a vaccine had been available after 100 days, over 71,000 deaths (16,507-204,249) and 384,000 (88,579-1,020,386) hospitalizations could have been averted. CONCLUSION: Our results highlight the substantial impact of NPIs, including lockdowns and curfews, in controlling the COVID-19 pandemic. We also demonstrate the value of the 100 days objective of the Coalition for Epidemic Preparedness Innovations (CEPI) initiative for vaccine availability.

Canadian Covid-19 pandemic public health mitigation measures at the province level.

The Covid-19 pandemic has prompted governments across the world to enforce a range of public health interventions. We introduce the Covid-19 Policy Response Canadian tracker (CPRCT) database that tracks and records implemented public health measures in every province and territory in Canada. The implementations are recorded on a four-level ordinal scale (0-3) for three domains, (Schools, Work, and Other), capturing differences in degree of response. The data-set allows the exploration of the effects of public health mitigation on the spread of Covid-19, as well as provides a near-real-time record in an accessible format that is useful for a diverse range of modeling and research questions.

Seasonality as a driver of pH1N12009 influenza vaccination campaign impact.

Although the most recent respiratory virus pandemic was triggered by a Coronavirus, sustained and elevated prevalence of highly pathogenic avian influenza viruses able to infect mammalian hosts highlight the continued threat of pandemics of influenza A virus (IAV) to global health. Retrospective analysis of pandemic outcomes, including comparative investigation of intervention efficacy in different regions, provide important contributions to the evidence base for future pandemic planning. The swine-origin IAV pandemic of 2009 exhibited regional variation in onset, infection dynamics and annual infection attack rates (IARs). For example, the UK experienced three severe peaks of infection over two influenza seasons, whilst Australia experienced a single severe wave. We adopt a seasonally forced 2-subtype model for the transmission of pH1N12009 and seasonal H3N2 to examine the role vaccination campaigns may play in explaining differences in pandemic trajectories in temperate regions. Our model differentiates between the nature of vaccine- and infection-acquired immunity. In particular, we assume that immunity triggered by infection elicits heterologous cross-protection against viral shedding in addition to long-lasting neutralising antibody, whereas vaccination induces imperfect reduction in susceptibility. We employ an Approximate Bayesian Computation (ABC) framework to calibrate the model using data for pH1N12009 seroprevalence, relative subtype dominance, and annual IARs for Australia and the UK. Heterologous cross-protection substantially suppressed the pandemic IAR over the posterior, with the strength of protection against onward transmission inversely correlated with the initial reproduction number. We show that IAV pandemic timing relative to the usual seasonal influenza cycle influenced the size of the initial waves of pH1N12009 in temperate regions and the impact of vaccination campaigns.

Pharmaceutical and non-pharmaceutical interventions for controlling the COVID-19 pandemic.

Disease spread can be affected by pharmaceutical interventions (such as vaccination) and non-pharmaceutical interventions (such as physical distancing, mask-wearing and contact tracing). Understanding the relationship between disease dynamics and human behaviour is a significant factor to controlling infections. In this work, we propose a compartmental epidemiological model for studying how the infection dynamics of COVID-19 evolves for people with different levels of social distancing, natural immunity and vaccine-induced immunity. Our model recreates the transmission dynamics of COVID-19 in Ontario up to December 2021. Our results indicate that people change their behaviour based on the disease dynamics and mitigation measures. Specifically, they adopt more protective behaviour when mandated social distancing measures are in effect, typically concurrent with a high number of infections. They reduce protective behaviour when vaccination coverage is high or when mandated contact reduction measures are relaxed, typically concurrent with a reduction of infections. We demonstrate that waning of infection and vaccine-induced immunity are important for reproducing disease transmission in autumn 2021.

Immunogenicity of COVID-19 vaccines and their effect on HIV reservoir in older people with HIV.

Older individuals and people with HIV (PWH) were prioritized for COVID-19 vaccination, yet comprehensive studies of the immunogenicity of these vaccines and their effects on HIV reservoirs are not available. Our study on 68 PWH and 23 HIV-negative participants aged 55 and older post-three vaccine doses showed equally strong anti-spike IgG responses in serum and saliva through week 48 from baseline, while PWH salivary IgA responses were low. PWH had diminished live-virus neutralization responses after two vaccine doses, which were 'rescued' post-booster. Spike-specific T cell immunity was enhanced in PWH with normal CD4+ T cell count, suggesting Th1 imprinting. The frequency of detectable HIV viremia increased post-vaccination, but vaccines did not affect the size of the HIV reservoir in most PWH, except those with low-level viremia. Thus, older PWH require three doses of COVID-19 vaccine for maximum protection, while individuals with unsuppressed viremia should be monitored for adverse reactions from HIV reservoirs.

Protocol for a living evidence synthesis on variants of concern and COVID-19 vaccine effectiveness.

BACKGROUND: It is evident that COVID-19 will remain a public health concern in the coming years, largely driven by variants of concern (VOC). It is critical to continuously monitor vaccine effectiveness as new variants emerge and new vaccines and/or boosters are developed. Systematic surveillance of the scientific evidence base is necessary to inform public health action and identify key uncertainties. Evidence syntheses may also be used to populate models to fill in research gaps and help to prepare for future public health crises. This protocol outlines the rationale and methods for a living evidence synthesis of the effectiveness of COVID-19 vaccines in reducing the morbidity and mortality associated with, and transmission of, VOC of SARS-CoV-2. METHODS: Living evidence syntheses of vaccine effectiveness will be carried out over one year for (1) a range of potential outcomes in the index individual associated with VOC (pathogenesis); and (2) transmission of VOC. The literature search will be conducted up to May 2023. Observational and database-linkage primary studies will be included, as well as RCTs. Information sources include electronic databases (MEDLINE; Embase; Cochrane, L*OVE; the CNKI and Wangfang platforms), pre-print servers (medRxiv, BiorXiv), and online repositories of grey literature. Title and abstract and full-text screening will be performed by two reviewers using a liberal accelerated method. Data extraction and risk of bias assessment will be completed by one reviewer with verification of the assessment by a second reviewer. Results from included studies will be pooled via random effects meta-analysis when appropriate, or otherwise summarized narratively. DISCUSSION: Evidence generated from our living evidence synthesis will be used to inform policy making, modelling, and prioritization of future research on the effectiveness of COVID-19 vaccines against VOC.

Immunogenicity of COVID-19 vaccines and their effect on the HIV reservoir in older people with HIV.

Older individuals and people with HIV (PWH) were prioritized for COVID-19 vaccination, yet comprehensive studies of the immunogenicity of these vaccines and their effects on HIV reservoirs are not available. We followed 68 PWH aged 55 and older and 23 age-matched HIV-negative individuals for 48 weeks from the first vaccine dose, after the total of three doses. All PWH were on antiretroviral therapy (cART) and had different immune status, including immune responders (IR), immune non-responders (INR), and PWH with low-level viremia (LLV). We measured total and neutralizing Ab responses to SARS-CoV-2 spike and RBD in sera, total anti-spike Abs in saliva, frequency of anti-RBD/NTD B cells, changes in frequency of anti-spike, HIV gag/nef-specific T cells, and HIV reservoirs in peripheral CD4 + T cells. The resulting datasets were used to create a mathematical model for within-host immunization. Various regimens of BNT162b2, mRNA-1273, and ChAdOx1 vaccines elicited equally strong anti-spike IgG responses in PWH and HIV - participants in serum and saliva at all timepoints. These responses had similar kinetics in both cohorts and peaked at 4 weeks post-booster (third dose), while half-lives of plasma IgG also dramatically increased post-booster in both groups. Salivary spike IgA responses were low, especially in INRs. PWH had diminished live virus neutralizing titers after two vaccine doses which were 'rescued' after a booster. Anti-spike T cell immunity was enhanced in IRs even in comparison to HIV - participants, suggesting Th1 imprinting from HIV, while in INRs it was the lowest. Increased frequency of viral 'blips' in PWH were seen post-vaccination, but vaccines did not affect the size of the intact HIV reservoir in CD4 + T cells in most PWH, except in LLVs. Thus, older PWH require three doses of COVID-19 vaccine to maximize neutralizing responses against SARS-CoV-2, although vaccines may increase HIV reservoirs in PWH with persistent viremia.

A pair formation model with recovery: Application to mpox.

The current global outbreaks of mpox is a unique infectious disease in the way it seems to be transmitting: it has been observed to be highly concentrated in communities of men who have sex with men (MSM) through pair formation, and also provides long lasting immunity. This framework of mostly close, prolonged contact spreading a disease that admits immunity after infection is unlike similar infections which either offer little to no immunity post-infection or are lifelong infections. This creates the need for a new model framework that incorporates pair formation structure with recovery. While seemingly a straight forward model, we show how new dynamics arise from the combination of pair formation and recovery that are not present in a standard model with recovery and also not present in a pair formation model without recovery. We see that the combination of these two properties allows for waves of infection that are not seen in a standard SIR model. These dynamics suggest that outbreaks of mpox around the world may require special attention from public health. We also derive a reproduction number for this model and estimate the reproduction number of human mpox to be ≈2.3 using global and Canadian data. The expression derived for R0 can help estimate key parameters for diseases transmission and public health interventions and compare to equivalent models without pair formation.

Monkeypox: a review of epidemiological modelling studies and how modelling has led to mechanistic insight.

Human monkeypox (mpox) virus is a viral zoonosis that belongs to the Orthopoxvirus genus of the Poxviridae family, which presents with similar symptoms as those seen in human smallpox patients. Mpox is an increasing concern globally, with over 80,000 cases in non-endemic countries as of December 2022. In this review, we provide a brief history and ecology of mpox, its basic virology, and the key differences in mpox viral fitness traits before and after 2022. We summarize and critique current knowledge from epidemiological mathematical models, within-host models, and between-host transmission models using the One Health approach, where we distinguish between models that focus on immunity from vaccination, geography, climatic variables, as well as animal models. We report various epidemiological parameters, such as the reproduction number, R0, in a condensed format to facilitate comparison between studies. We focus on how mathematical modelling studies have led to novel mechanistic insight into mpox transmission and pathogenesis. As mpox is predicted to lead to further infection peaks in many historically non-endemic countries, mathematical modelling studies of mpox can provide rapid actionable insights into viral dynamics to guide public health measures and mitigation strategies.

Determination of significant immunological timescales from mRNA-LNP-based vaccines in humans.

A compartment model for an in-host liquid nanoparticle delivered mRNA vaccine is presented. Through non-dimensionalisation, five timescales are identified that dictate the lifetime of the vaccine in-host: decay of interferon gamma, antibody priming, autocatalytic growth, antibody peak and decay, and interleukin cessation. Through asymptotic analysis we are able to obtain semi-analytical solutions in each of the time regimes which allows us to predict maximal concentrations and better understand parameter dependence in the model. We compare our model to 22 data sets for the BNT162b2 and mRNA-1273 mRNA vaccines demonstrating good agreement. Using our analysis, we estimate the values for each of the five timescales in each data set and predict maximal concentrations of plasma B-cells, antibody, and interleukin. Through our comparison, we do not observe any discernible differences between vaccine candidates and sex. However, we do identify an age dependence, specifically that vaccine activation takes longer and that peak antibody occurs sooner in patients aged 55 and greater.

Multiple cohort study of hospitalized SARS-CoV-2 in-host infection dynamics: Parameter estimates, identifiability, sensitivity and the eclipse phase profile.

Within-host SARS-CoV-2 modelling studies have been published throughout the COVID-19 pandemic. These studies contain highly variable numbers of individuals and capture varying timescales of pathogen dynamics; some studies capture the time of disease onset, the peak viral load and subsequent heterogeneity in clearance dynamics across individuals, while others capture late-time post-peak dynamics. In this study, we curate multiple previously published SARS-CoV-2 viral load data sets, fit these data with a consistent modelling approach, and estimate the variability of in-host parameters including the basic reproduction number, R0, as well as the best-fit eclipse phase profile. We find that fitted dynamics can be highly variable across data sets, and highly variable within data sets, particularly when key components of the dynamic trajectories (e.g. peak viral load) are not represented in the data. Further, we investigated the role of the eclipse phase time distribution in fitting SARS-CoV-2 viral load data. By varying the shape parameter of an Erlang distribution, we demonstrate that models with either no eclipse phase, or with an exponentially-distributed eclipse phase, offer significantly worse fits to these data, whereas models with less dispersion around the mean eclipse time (shape parameter two or more) offered the best fits to the available data across all data sets used in this work. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

COVID-19 Vaccination and Healthcare Demand.

One of the driving concerns during any epidemic is the strain on the healthcare system. As we have seen many times over the globe with the COVID-19 pandemic, hospitals and ICUs can quickly become overwhelmed by cases. While strict periods of public health mitigation have certainly helped decrease incidence and thus healthcare demand, vaccination is the only clear long-term solution. In this paper, we develop a two-module model to forecast the effects of relaxation of non-pharmaceutical intervention and vaccine uptake on daily incidence, and the cascade effects on healthcare demand. The first module is a simple epidemiological model which incorporates non-pharmaceutical intervention, the relaxation of such measures and vaccination campaigns to predict caseloads into the Fall of 2021. This module is then fed into a healthcare module which can forecast the number of doctor visits, the number of occupied hospital beds, number of occupied ICU beds and any excess demand of these. From this module, we can also estimate the length of stay of individuals in ICU. For model verification and forecasting, we use the four most populous Canadian provinces as a case study.

A mathematical model of protein subunits COVID-19 vaccines.

We consider a general mathematical model for protein subunit vaccine with a focus on the MF59-adjuvanted spike glycoprotein-clamp vaccine for SARS-CoV-2, and use the model to study immunological outcomes in the humoral and cell-mediated arms of the immune response from vaccination. The mathematical model is fit to vaccine clinical trial data. We elucidate the role of Interferon-γ and Interleukin-4 in stimulating the immune response of the host. Model results, and results from a sensitivity analysis, show that a balance between the TH1 and TH2 arms of the immune response is struck, with the TH1 response being dominant. The model predicts that two-doses of the vaccine at 28 days apart will result in approximately 85% humoral immunity loss relative to peak immunity approximately 6 months post dose 1.

The race to understand immunopathology in COVID-19: Perspectives on the impact of quantitative approaches to understand within-host interactions.

The COVID-19 pandemic has revealed the need for the increased integration of modelling and data analysis to public health, experimental, and clinical studies. Throughout the first two years of the pandemic, there has been a concerted effort to improve our understanding of the within-host immune response to the SARS-CoV-2 virus to provide better predictions of COVID-19 severity, treatment and vaccine development questions, and insights into viral evolution and the impacts of variants on immunopathology. Here we provide perspectives on what has been accomplished using quantitative methods, including predictive modelling, population genetics, machine learning, and dimensionality reduction techniques, in the first 26 months of the COVID-19 pandemic approaches, and where we go from here to improve our responses to this and future pandemics.

An HIV stochastic model with cell-to-cell infection, B-cell immune response and distributed delay.

In this study, a delayed HIV stochastic model with virus-to-cell infection, cell-to-cell transmission and B-cell immune response is proposed. We first transform the stochastic differential equation with distributed delay into a high-dimensional degenerate stochastic differential equation, and then theoretically analyze the dynamic behaviour of the degenerate model. The unique global solution of the model is given by rigorous analysis. By formulating suitable Lyapunov functions, the existence of the stationary Markov process is obtained if the stochastic B-cell-activated reproduction number is greater than one. We also use the law of large numbers theorem and the spectral radius analysis method to deduce that the virus can be cleared if the stochastic B-cell-inactivated reproduction number is less than one. Through uncertainty and sensitivity analysis, we obtain key parameters that determine the value of the stochastic B-cell-activated reproduction number. Numerically, we examine that low level noise can maintain the number of the virus and B-cell populations at a certain range, while high level noise is helpful for the elimination of the virus. Furthermore, the effect of the cell-to-cell infection on model behaviour, and the influence of the key parameters on the size of the stochastic B-cell-activated reproduction number are also investigated.

Assessing transmission risks and control strategy for monkeypox as an emerging zoonosis in a metropolitan area.

To model the spread of monkeypox (MPX) in a metropolitan area for assessing the risk of possible outbreaks, and identifying essential public health measures to contain the virus spread. The animal reservoir is the key element in the modeling of zoonotic disease. Using a One Health approach, we model the spread of the MPX virus in humans considering potential animal hosts such as rodents (e.g., rats, mice, squirrels, chipmunks, etc.) and emphasize their role and transmission of the virus in a high-risk group, including gay and bisexual men-who-have-sex-with-men (gbMSM). From model and sensitivity analysis, we identify key public health factors and present scenarios under different transmission assumptions. We find that the MPX virus may spill over from gbMSM high-risk groups to broader populations if the efficiency of transmission increases in the higher-risk group. However, the risk of outbreak can be greatly reduced if at least 65% of symptomatic cases can be isolated and their contacts traced and quarantined. In addition, infections in an animal reservoir will exacerbate MPX transmission risk in the human population. Regions or communities with a higher proportion of gbMSM individuals need greater public health attention. Tracing and quarantine (or "effective quarantine" by postexposure vaccination) of contacts with MPX cases in high-risk groups would have a significant effect on controlling the spreading. Also, monitoring for animal infections would be prudent.